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Abstract 

The replacement paths problem for directed graphs is to find for given nodes s and t and every edge 
e on the shortest path between them, the shortest path between s and t which avoids e. For unweighted 
directed graphs on n vertices, the best known algorithm runtime was 0(n 2 ' 5 ) by Roditty and Zwick. For 
graphs with integer weights in {— M, . . . , M}, Weimann and Yuster recently showed that one can use fast 
matrix multiplication and solve the problem in 0(Mn ' ) time, a runtime which would be 0(Mn ' ) 
if the exponent ui of matrix multiplication is 2. 

We improve both of these algorithms. Our new algorithm also relies on fast matrix multiplication 
and runs in 0(Mn^ poly log n) time if w > 2 and 0(Mn +e ) for any e > if lo — 2. Our result shows 
that, at least for small integer weights, the replacement paths problem in directed graphs may be easier 
than the related all pairs shortest paths problem in directed graphs, as the current best runtime for the 
I] , latter is 0(n 2,5 ) time even if uj = 2. 
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1 Introduction 

In the replacement paths problem, one is given a graph and a shortest path P between two vertices s and 
^O ' t, and one needs to return for every edge on P, the shortest path between s and t which avoids e. The 

typical motivation for the problem is that links in a network can fail and a backup path between important 
vertices would be useful. Replacement paths are useful also in contexts in which one may wish to satisfy 
other constraints beyond short length |15j . For instance, in biological sequence alignment |2J replacement 
paths are useful in determining which pieces of an alignment are most important. The replacement paths 
problem is also used in the computation of Vickrey Prices of edges that are owned by selfish agents in a 
network [15] Q2 , and in finding the k shortest simple paths between two nodes [TSl E£H HH1 120] • 

The replacement paths problem for undirected graphs can be solved very efficiently: Malik et al. [TB] gave 
an 0(mj^| algorithm for graphs on m edges; Nardelli et al. [17] used Thorup's linear time algorithm for single 
source shortest paths to improve the runtime to 0(ma(n)) in the word- RAM model of computation. The best 
algorithm for the problem in directed graphs with arbitrary edge weights is by Gotthilf and Lewenstein |llj 
and runs in 0(mn + n 2 log log n) time, where m and n are the number of edges and vertices respectively. 
For dense graphs, nothing much better than cubic time is knowno Recently, Vassilevska- Williams and 
Williams [22] showed that the general replacement paths problem in directed graphs is equivalent to the 
all pairs shortest paths problem (APSP), under subcubic reductions. This essentially means that either 
both problems admit truly subcubic algorithms, or neither of them does. It is worth pointing out that this 
apparent cubic time barrier is only due to the wish to compute the replacement distances exactly. In contrast, 
Bernstein [2] recently showed that there is an almost linear time approximation scheme for replacement paths. 
For unweighted directed graphs, Roditty and Zwick 20 gave a randomized combinatorial algorithm which 
computes replacement paths in 0(m^/n) time. In their recent FOCS paper, Weimann and Yuster [21] were 
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the first to apply fast matrix multiplication techniques to the problem. Using the current best algorithm 
for matrix multiplication [6], Weimann and Yuster obtain a randomized algorithm for replacement paths 
which runs in 0(M n 2 ' 584 ) for n node graphs with integer weights in {— M, . . . , M}. It has the advantage 
that it can also be used to construct a distance sensitivity oracle. Furthermore, if the exponent u> of matrix 
multiplication is actually 2, their runtime would be 0(Mn 233 ) which would improve on Roditty and Zwick's 
0(n 25 ) runtime for unweighted graphs. 

In this work we reconsider the replacement paths problem for dense directed graphs with weights in 
{— M, . . . , M}. We give a simple reduction from the problem to all pairs shortest paths (APSP) and by 
applying Zwick's algorithm [26] , we obtain an improved deterministic algorithm for replacement paths. 

Theorem 1. There is a deterministic algorithm for the replacement paths problem in directed graphs with 
weights in {— M, . . . , M } which runs in 0(M°' 681 n 2 ' 575 ) time. 

This runtime bound improves on Weimann and Yuster's randomized algorithm (although if u = 2, 
Weimann and Yuster's bound would be lower). We then combine the ideas from the APSP reduction with 
a divide and conquer approach to design a faster randomized algorithm for replacement paths. 

Theorem 2. The replacement paths problem in a directed graph with integer weights in {— M, . . . , M} can 
be computed by a randomized algorithm in O(M ' 832 n 2 ' 481 ) time. 

Finally, we use a recursive approach on top of the ideas used in proving Theorem[5]to improve the runtime 
even further: 

Theorem 3. Let w < 2.376 JS^ be the exponent for matrix multiplication. The replacement paths problem 
in a directed graph with integer weights in {— M, ...,M} can be computed by a randomized algorithm in 
O ( Mn u poly log n) time, if u> > 2, or in 0(Mn 2+e ) for every e > if to = 2. 

Because of prior reductions between the replacement paths problem and the k shortest simple paths 
problem [50] we obtain the following corollary. 

Corollary 1. The k shortest simple paths in a directed graph with edge weigths in {—M,...,M} can be 
computed in randomized 0(k ■ Mn"poly log n) time. 

Theorem [3] improves on Roditty and Zwick's 0(n 25 ) runtime. It also improves the range of M for which 
there is a subcubic algorithm for the problem: the previous range of Weimann and Yuster's algorithm was 
for M — O(n 0Ale ), and our algorithms are subcubic for any M — O(n 0623 ). Our result also shows that, at 
least for small integer weights, the replacement paths problem in directed graphs might actually be easier 
than APSP in directed graphs. The current best algorithm for APSP in unweighted directed graphs is by 
Zwick [26] : it uses fast matrix multiplication and would run in J7(n 25 ) time even if w = 2. Furthermore, 
improving on 0(n w ) for replacement paths in directed unweighted graphs would likely require radically new 
techniques, as the problem is very related to Boolean matrix multiplication. 

Other related work. The replacement paths problem is closely related to the problem of finding the 
second shortest path between two nodes, and in general to the k shortest paths problem. For directed graphs 
with m edges and n vertices and nonnegative edge weights, Eppstein [8] gave an algorithm which returns 
the k shortest paths between two given nodes s and t in time 0(m + nlogn + k) . The paths that Eppstein's 
algorithm returns, however, may not be simple. When the k shortest paths are required to be simple, the 
fastest known algorithms for the problem use algorithms for the replacement paths problem. Roditty and 
Zwick [50] showed that the k simple shortest paths problem can be reduced to 0(k) computations of the 
second shortest simple path. For unweighted graphs, both the replacement paths problem and the second 
shortest simple path problem are closely related to Boolean matrix multiplication; for general directed 
graphs with arbitrary weights, both problems are equivalent to all pairs shortest paths under subcubic 
reductions [22]. Hershberger et al. 13 showed that for directed graphs with arbitrary edge weights, the 
replacement paths problem requires J7(n 25 ) time in the path-comparison model of computation of Karger 
et al. IH. 



Preliminaries. All graphs in this paper have n vertices, unless otherwise stateds. They are directed and 
weakly connected. Each graph G — (V, E) is equipped with a weight function, £ : E — > {—M, . . . , M} where 
M is a positive integer. The length or weight of a path v± — > i>2 —>...—> vt in G is just X)i=i ^( u jj w «+i)- A 
shortest path between two nodes is just a path between them which minimizes the length function. 

Given a path P = vi —¥ v-i — > . . . — > v t in G, let dp(vi,Vj) = ^ZiZ^ £(vk,Vk+i)^ i-6. the length of the 
subpath of P from Uj to Uj . 

For an integer n > 0, let [n] = {1, 2, . . . , n}. Let A be an m x n matrix and B be an n x p matrix, where 
A and i? have entries from ZU {oo}. Then their distance product A* B is defined as 

(A*P)[i,j] = mm(A[i,k]+B[k,j}), V^ € [m], j € [p]. 
fce[n] 

Finally, let w be the infimum of all numbers such that there exists a constant uq so that for all n > no, 
the product of two n x n integer matrices can be computed in 0{n u ) time. 

Replacement paths problem. Given G = (V, E) with edge weights in {— M, . . . , M} and s,t £ V, the 
task is to compute the shortest path P from s to t in G and for every edge e £ P, the shortest path in 
G \ {e}. It is known how to compute the actual shortest path P efficiently For instance, when the weights 
are all nonnegative, one can just use Dijkstra's algorithm [7] with Fibonacci heaps [9] and find P in 0(m + 
n log n) time which is linear for dense graphs. When the weights can be negative, Goldberg's algorithm [TO] 
shows how to find P or detect a negative cycle in the graph in 0(m^/n poly log M) time. Yuster and 
Zwick's [24] algorithm relies on fast matrix multiplication and finds P in 0(Mn u ') time; using Coppersmith 
and Winograd's [5] bound on u, their running time is 0(Mn 2376 ) which is better than Goldberg's runtime 
for dense graphs and small M. Relying on the existing subcubic algorithms, we can from now on assume 
that P is given as part of the input. 

The naive way to compute the replacement paths is to remove each edge e £ P in turn and compute the 
shortest path in G \ {e} from scratch. This approach is however unnecessarily time intensive - the shortest 
paths computations share a lot of information. 

2 Two Simple Algorithms 

The notion of a detour is very important for the replacement paths problem. Let P — {s = v± — > vi — >...—> 
Vk = i] be the shortest path between s and t in G. Let E{P) denote the set {(uj,Uj+i)}ig[fc_ 11 of all edges 
on P. For j > i, a detour D(vi,Vj) between Vi and Vj is a simple path from Vi to Vj which does not contain 
any other nodes of P. A detour D(vj, Vk) is said to circumvent edge (vi, Vj+i) if j < i and i + I < k. It is 
well known (see e.g. [2], Lemma 2.1) that for any edge e = (WjjUi+i) £ E(P), the shortest path between s 
and t in G \ {e} is exactly the minimum out of all paths of the form 

s ->• v% -¥ . . . -*• Vj D(vj , Vk) Vk -> . . . ->• t, 

where j < i and i + 1 < k. 

From the above discussion, we see that in order to find all replacement paths, it is sufficient to compute 
all the detours between nodes in P. After all detours are found, we can compute in 0(n 2 ) extra time for 
every detour D(vi,Vj), the length £(D(vi, Vj)) of the path a —>■ . . . —> Vj D(vj, Vk) Vk —> ■ ■ ■ — >• t which it 
defines. This can be done by first computing the weight of each subpath s — >• . . . — > Vi for all i and the weight 
of each subpath Vj —>...—>• t, via dynamic programming; then £(D(vi,Vj)) can be computed in constant 
time by adding the weights of the two paths and that of the detour. 

After this we can sort the detours in nondecreasing order according to £(■) in 0(n 2 logn) time. We can 
then store all nodes Vi of P in a successor search data structure T (e.g. any balanced binary search tree), 
sorted according to i, their position in P. To process the current shortest detour D(vi,Vj) in the sorted 
order, find the successor v x of Vi in T, and if a; < j — 1, record that the shortest replacement path length_| 
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Figure 1: An example of the reduction from replacement paths to APSP in the case of an unweighted graph. 

for (v x ,v x -\-\) is £(D(vi,Vj)). Then remove v x from T and process the rest of the elements in T which are 
between Vi and Vj. Collecting all the results from the detour information only takes an additional 0(n 2 logn) 
time. The algorithms below show how to compute the detours. 

Algorithm 1: A reduction to APSP. We show that the detours can be computed using a reduction to 
APSP. The possibility of such a reduction was mentioned in the work of Bhosle and Gonzalez [3J; here we 
make it explicit and show that it can be used to further improve the known runtime bounds. 

For every node Vi of P, create two copies v™ and v° ut . Create a new graph G' by taking G \ P U 



K 



i v i >ie[fc]- 



For every edge (vi,u) for u £ G \ P, add an edge (u| 



of the same weight in G' . 



Similarly, for every edge (u, Vi) for u E G\P, add an edge (it, v\ n ) of the same weight in G' . G' is essentially 
G with the edges of P removed, except that each node of P is split into two. See Figure [TJ for an example 
conversion from G to G'. Now compute all pairs shortest paths in G". The shortest path between v° ut and 
v 1 ™ is exactly the optimal detour D(vi, Vj) in G. Thus with one call to an APSP algorithm, and 0(n 2 logn) 
extra time we obtain an algorithm for replacement paths. 

Theorem 4 (Zwick's APSP |26j). The APSP problem in a directed graph with integer edge weights in 
{—M, . . . , M} can be computed in 0(At 1 ' ( 4_LJ )n 2+1 ' ( 4 ~ w ') time, where u> < 2.376 is the exponent of matrix 
multiplication. If rectangular matrix multiplication is used, then the running time is O(A/ 0681 n 2575 ). 

Applying Zwick's [26] algorithm for APSP from Theorem [4] with weights in {— M, ...,M} we prove 
Theorem [TJ 

Reminder of Theorem [TJ There is a deterministic algorithm for the replacement paths problem in 
directed graphs with weights in {—A/, . . . , M} which runs in O(M ' 681 n 2 ' 575 ) time. 



Algorithm 2: Divide and Conquer. This second algorithm has a worse running time, but illustrates 
the fact that the problem lends itself to some partitioning. Consider a subpath P' of P lying between nodes 
v x and v y . Let V(P') be the nodes of P' except for v x and v v . Consider an edge (v^, Wi+i) in P'. Then, the 
best detour D(vj,Vk) which circumvents (vi,Vi+i) either has both Vj and Vk in V(P'), or Vj is to the left of 



V(P') and v k is in V(P') 



is in V(P') and v k is to the right ot V(P'), or Vj is to the left and v k is to 



the right of V(P') and D(vj,Vk) does not touch V(P') at all. 

Suppose we have precomputed in 0(n 2 ) time, for every Vi and Vj the length dp(vi, Vj) of the subpath of 
P from Vi to Vj. Consider any subpath P' as above. We can create a graph Gp* as follows. Let P l and P r 
be the subpaths of P to the left and right of P' respectively. Take G and remove all incoming edges to nodes 
on P l except those in P and all outgoing edges from nodes on P r except those in P. This will ensure that 
any path that we compute exiting P l doesn't reenter it and any path entering P r doesn't rccxit it. Now 
remove all edges in P' and split each node v in V(P') into two as before: a copy v m with all the remaining 
incoming edges and a copy v out with all the remaining outgoing edges. An example of the construction of 
Gp> is given in Figure [5] 
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Figure 2: An example of the transformation of a graph G into Gp>, given P and subpath P'. Notice that 
edges (y, 1*2) and (w5,y) are removed, and that all nodes in V(P') get split into two. The relevant paths in 
Gp' corresponding to detour paths circumventing edges in P are s — > V2 — > X — > v™, s ^- V2 ^- x — > y — > 
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Consider any (simple) path from v° ut to vj 1 for j > i and Vi,Vj G V(P'). It corresponds to a detour in 
G since it does not use any other nodes in P by construction. Consider any path from if 11 * £ V(P') to t in 
Gp> . It corresponds to a portion of some replacement path for the edges of P' after Vi . If we know the length 
of the shortest such path between v° ut and t, then we can compute the length of the shortest replacement 
path in G that follows P from s to Vi and circumvents all edges in P' after Vi . Similarly, if we know the 
length of the shortest such path between s and uj" in Gpi , then we can compute the length of the shortest 
replacement path in G that circumvents all edges in P' before Vi and follows P from Vi to t. Finally, the 
length of the shortest path between s and t in Gpi is exactly the length of the shortest path in G which 
circumvents all edges in P'. Thus, if we can compute the shortest paths in Gpi between all pairs of nodes 
among the copies of nodes in V(P'), s and t, then we can find the replacement paths for all edges in P' in 
0(|P'| 2 logn) extra time with the successor search data structure approach from before. To compute these 
paths for P' we use the following theorem by Yuster and Zwick [53] . 

Theorem 5 (Yuster-Zwick |24j). Given a directed graph G with edge weights in [— M, ...,M], one can 
compute in 0{Mn u) ) time an n x n matrix D, so that the (i,j) entry of the distance product D * D is the 
distance between nodes i and j in G. 

Since we only need to compute 0(|P'| 2 ) pairs of distances in Gpi, given the matrix D for Gpi from 
Theorem [5j we can do this in overall 0{Mn u + n\P'\ 2 ) time. If \P'\ = v'ln ( ^ 1,/2 , then the runtime is 
6{Mn"). _ 

This suggests the following bucketting approach: partition P into q < n^ 3_ "^ 2 /-\/M pieces Pi, . . . ,P q 
of size roughly n/q — -s/Mn^ -1 " 2 . Create a graph Gp t for each piece i £ [q] as before. Compute the 
relevant shortest paths in overall time 0(qMn u ) — 0(VMn^ 3+ul ^ 2 ) . After this, compute the lengths of all 
viable replacement paths in 0(n 2 ) extra time, and finally find the best replacement path for each edge using 
successor data structures in 0(n 2 log n) time. The overall runtime of the algorithm is 0(vMn' 3+IJ " 2 ) = 
0(\/Mn 2 ' 6S8 ), slightly slower than the APSP-reduction technique in Algorithm 1. 



3 A faster algorithm: Algorithm 3. 

Here we will combine the above two algorithms to prove Theorem [5] The idea is as follows: we use some 
random sampling and Theorem [5] combined with fast matrix multiplication to compute partial distances 
between roughly equally spaced nodes on all detours. After this the bucketting technique of Algorithm 2 
can actually be used to reduce the replacement paths problem to n/b instances of APSP on graphs with 
0(b) nodes and with only slightly increased edge weights. We can then use Zwick's APSP algorithm on 
the pieces and afterwards combine the results in roughly quadratic time. In the following we will ignore 
polylogarithmic factors in the runtime. 



Partitioning all long detours. Call any detour on at least n 1 ~ 6 edges long. Our algorithm will be able 
to handle short detours easily. Here we split the long detours into pieces which we will be able to handle just 
like the short detours. Following an idea from [2T] we note that if we pick a set B of size roughly b — Cn 6 In n 
uniformly at random for large enough constant C, then every long detour path D is partitioned by points of 
B into pieces of size at most n 1 ^ 6 with high probability. 

Lemma 1. Let N > n. Let B be a random sample of (C + 2>)n e h\N nodes ofV. With probability at least 
1 — 1/N C , the nodes of B touch all detours on at least n 1_e edges, so that the nodes of B partition every 
detour into subpaths on at most 2n 1 ~ £ edges each. 

Proof. Let Q be a collection of n 3 paths, each on n 1 ^ 6 nodes. We first show that the probability that B 
contains at least one node from each path in Q is at least 1 — 1/N C . We can assume that B was formed by 
picking each node independently at random with probability 1/n. Let P £ Q. The probability that none of 
the nodes of B are in P is (1 - l/n £ )( c+3 )" e ln N < 1/N C+3 . The probability that there exists a path in Q 
such that none of its nodes are in B is 1/N by a union bound. 

Now consider each detour v\ —>...—> v t on at least n 1_e nodes. Partition it by splitting off subpaths on 
n x ~ e edges from one end until there is one path on < n l ~ E nodes left. This creates 0(n e ) disjoint subpaths 
on n 1_£ nodes each and possibly one short path. Let Q be the set of all long subpaths of all 0(n 2 ) long 
detours. There are at most n 3 of these and by our above argument, every one of them is hit by B with high 
probability. Hence the number of edges between consecutive points of B in each detour is at most 2n 1 ~ e , 
and for any detour the number of edges from its beginning and the first point of B on it and from its end 
and the last point of B is also at most 2n 1 ~ e . □ 

Computing short detour subpaths. We compute the matrix D from Theorem [3] for the graph G' used 
in Algorithm 1 (where the edges of P are removed and the nodes are split into two). We only care about 
the distances between nodes of B and nodes of B U P in G' which are of long detours, and about distances 
between nodes of P which use short detours. All such distances are less than 0(Mn 1 ~ e ). Hence, we can 
disregard any entries of D greater than 0(Mn}~ e ). Then we can use the following theorem of Alon, Galil 
and Margalit p] (following Yuval [25]). 

Theorem 6 (1, 25]). One can compute the distance product of two nxn matrices with entries in {~N, . . . , N} 
in 6(Nn u ) =0(Nn 23re ) time. 

Using Theorem HI we can compute the lengths w(-) of all short detour paths and all paths between nodes 
of B and nodes of B U P lying on long detours. The runtime is 0(Mn 1 ^ e+u ) = 0(Mn 3 - 376 ~ £ ). 

Pasting the subpaths together to compute replacement paths. We set aside the short detours. We 
process them in additional 0(n 2 logn) time just as in Algorithm 1. Now we handle the rest of the detours. 
Since we know that the nodes of B partition the long detours and we have computed all short paths between 
nodes of B and nodes of B U P, we can focus on the graph which is only composed ot B U P with the 
new weights w(-) corresponding to short (on < 2n 1 ~ e edges) pieces of long detours. All the weights are 
0(Mn l ~ e ). Now we use the partitioning idea of Algorithm 2 to split the problem into 0{n/\B\) pieces each 
of size 0(|B|) = 0(n e ). If we use Algorithm 2 naively, even though B is small, we would actually be left 
with a large part of P in each subinstance, and the subinstance size would be linear in n. This can be fixed, 
however. 

As in Algorithm 2, we partition P into pieces {Pi, . . . , P n /\B\} of si ze roughly \B\. We then create graphs 
Gp t for each piece. Consider a subinstance Gp i . The nodes corresponding to nodes of P are just the subpaths 
P l and P r to the left and right of P, with edges into P l and out of P r removed, together with in and out 
copies of the nodes of Pj. The number of nodes corresponding to Pi is already 0(|P|). We will argue that 
we can replace P l and P r by 0(|P|) nodes each so that the edge weights of the graph increase by at most 
0(Mn ) and the distances between the relevant node pairs are preserved. 

We will show how to replace P r ; P l can be processed similarly Let P r = {pi, . . . ,pt} for t < n. Split P r 
into t/c = 0(n E ) subpaths {{pi -» p c ), (p c +i -> P2c), • • • , (p((t/ c )-i)c+i -> Pt)} on c - 1 = n 1_e consecutive 





Figure 3: An example of the transformation of path P r in our main algorithm. Each consecutive subpath 
gets a node Ui corresponding to the last node on the subpath. The nodes Ui are chained in a path with 
weights corresponding to the distance of the last subpath nodes in the original path P r . This is possible 
since any path that we are interested in never leaves P r once it enters it. 



nodes. We will replace for each i > 0, the subpath (pi C +i — > Pu+i) c ) °y a new node u^+i so that for every 
node b 6 B, the edge weight of (6,-Ui+i) is min°~ w(b,pi C +i+j) + dp(pi C +i+j , P(i+i) c ) ■ As in Algorithm 1, 
there are no edges from any Ui to B. Connect the Ui nodes by a path by adding edges (ui, Wi+i) with weight 
dp(j>i C -,P(i+i)c) f° r all 1 < i < t/c. See Figure [3] for an example transformation of a path P r . 

The weights of the edges in the new graph are still 0(Mn 1 ~ e ). The distance between any node in Pi and 
t is the same as before. The distance between s and t is the same as before. Assuming that all distances 
dp(x,y) are precomputed, the time to do this transformation is 0(n|_B|) per piece Pi, and is thus 0(n 2 ) 
overall. 

Thus, we have partitioned the problem in 0(n 2 + Mn 3,3 ™~ £ ) time into 0{n/\B\) = 0(n l ~ E ) instances of 
APSP on graphs with 0(\B\) = 6{n e ) nodes with weights between {-N, . . . , N} where N = 0(Mn 1 - £ ). In 
each subinstance, we can compute APSP using Zwick's algorithm [26] in time 0((Ajfn 1-e )°- 681 (n £ ) 2 - 575 ) = 



0(M 



0.681^0.681 + 1. 894 E 



). The running time of this step over all 0(n £ ) subinstances is asymptotically 



'■[M 



, r 0. 681^0. 681^ 



1.894£i 



M 



0.681^1. 681+0. 894e 



The overall runtime is 



7Vf°-681 n l. 681+0. 894e 



Mn 3.376-e 



which is minimized when M - 681 ™ 1 - 68 ^ - 8946 = Mn 3 - 376 ^, 



when n 1894£ = M a319 n L695 . Then rf 



M°- 



168^0.895 



and the runtime is 0(M 



0.832^2.481 



). We have proven Theorem [2] 



4 An even faster algorithm. 

In the previous section we reduced in 0(Mn l+u ~ £ ) the replacement paths problem on n-node graphs with 
weights in {— M, . . . , M} to n 1 ^ 6 instances of APSP in graphs on n e nodes with weights in {— Mn 1 ^ 6 , . . . , Mn 1 ^ 6 }. 
Here we make the observation that the APSP computation is not necessary, and we can instead solve the 
problem using recursion. 

We will solve instances of the following problem. 

Definition 4.1 (Circumventing Paths Problem). Let the following be given: Z,n,p,M , and T = n/z p , a 
graph G on 0(T\og p n) nodes with edge weights in {—N, . . . ,N} with N = 0(MZ p ), an ordered list of T 
tuples ((u™, ii° ut ), . . . , (u|™, i>°"')), and for every j G [T], values dp(s,Vj) and dp(vj,t). Every v™ has only 
incoming edges and each v° ut has only outgoing edges. Compute for every x S [T] : 

min d P (s,v h ) + d(v°f,v%) + d P {v lk ,t), 



where d(v° 



, v]™) is the distance in G between v° ut and v] n . 



Each ordered list of tuples above will correspond to a subpath P' of P, and dp(u, v) is just the distance 
between u and v inside P. Now consider an instance of Circumventing Paths as in the definition and let P' 
be a path Vi 1 — >• . . . — > t>i T corresponding to the ordered list ((u™, f°"*), • ■ • , (u| n , v° ut )). Consider a subpath 
P, of P' on T/Z = n/Z p+1 nodes. Let P be a random sample of (J(T log n) nodes from G\P' so that all 
paths between pairs of nodes on Pi on at least Z nodes are partitioned by B into pieces of length O(Z), as 
in Lemma [I] Just as in Algorithm 3, use Theorems [3] and [5] and compute all distances shorter than MZ P+1 
between nodes of G. This finds all detours on a small number of nodes. We can then compute the length of 
the paths in the original graph G that these detours create and add the corresponding replacement paths to 
a global list L. 

After all this computation, we can create two graphs G(Pi) and G Pi for Pi similar to Algorithm 3. 

The graph G(Pi). The graph G{Pi) for Pj is created as follows. First we take P' U B and put an edge 
between any two nodes of B with a weight corresponding to their distance in G, if it was computed. We 
also add an edge between an out- node of Pi and a node of B, or a node of B and an in- node of Pi, using the 
distance (if it was computed) as a weight. The weights of the graph are now 0(M Z p+1 ). 

As before, let P l and P r denote the subpaths of P' to the left and right of Pi, respectively, where 
P' = P l Pi P r (0 stands for concatenation). Abusing notation somewhat, let P l and P r also correspond 
to the copies of nodes of P l and P r which are in the ordered list of tuples of the Circumventing Paths 
instance. 

We remove all in-copies of nodes in P l and all out-copics of nodes of P r . Then, just as in Algorithm 
3, we replace P l and P r each with a path on roughly T/Z nodes, corresponding to subpaths on roughly Z 
nodes each. We define new weights in these paths and between B and the nodes of these paths according to 
distance, as in Figure [3] The main difference is that Vi x plays the role of s and Vi T plays the role of t. More 
formally, we process P r as follows (P l is processed symmetrically): Let P r = {p\, . . . ,pt}, where p t — w, T . 
Split P r into 0{t/z) = 0{T/Z) subpaths {(p x -> p c ), (p c+1 -> p 2c ), ■■■, (p((t/o)-i)c+i ~> Pt)} on c - 1 = Z 
consecutive nodes. We will replace for each i > 0, the subpath (pi C +i — > Pu+i) c ) by a new node Uj+i so 
that for every node b £ P, the edge weight of (6, Uj+i) is min^~ w(b,pi C+ i + j) + dp(pi c +i+j,P(i+i)c), where 
w(-, •) are the distances computed in G. Connect the Ui nodes by a path by adding edges (ui,Ui+±) with 
weight dp(pic,p( i+ i\ c ) for all 1 < i < t/c. This increases the maximum edge weight in the graph by at most 
0{MZ p+1 ). 

In each graph G(Pi) we only need to compute: 

1. the shortest path from the node corresponding to Vi 1 to all in-copies of nodes of V(Pj); 

2. the shortest path from all out-copies of nodes of V(Pi) to the node corresponding to Vi T ; 

3. the shortest path between the node corresponding to Vi t to the node corresponding to Vi T ; 

4. for every r £ [T/Z], the length I of a path between v° ut and u|" in G(Pi) with j < r < k, such that 
the quantity dp(s, Vi j ) + £ + dp(i)i k , t) is minimized over j, k with j < r < k such that Vi j , Vi k £ Pi. 

1 and 2 compute the shortest pieces of detours one of the endpoints of which is not in V(Pj). 3 computes 
an optimal detour circumventing all edges in Pi, and 4 is concerned with detours with endpoints within 
V(P). 

We note that one can compute 1, 2 and 3 in time 0((M Z P+1 )(T \ogn/ Z)^) using Yuster and Zwick's single 
source shortest paths algorithm [23]. After those distances are computed, one can convert them to lengths of 
candidate replacement paths as follows: Suppose d is a distance computed between the node corresponding 
to Vi x and a node %)%, £ V(Pj). Then dp{s,Vi t ) + d + dp{vu ,t) is the length of some replacement path which 
circumvents all edges of P* before Vi-. Similarly, if d is the distance computed between a node Vi- £ V(Pj) 
and the node corresponding to u, T , then dp(s,Vi) + d + dp(vi T , t) is the length of some replacement path 
which circumvents all edges of Pi before Vi . . 



The graph G Pi . The graph G Pi is formed from G(Pt) by removing the nodes corresponding to P l and P r . 
Point 4 from above can be done with recursive calls to computing circumventing paths in the graph G Pi . 

Pseudocode for our approach is given in Algorithm [1] and Procedures [5] and [3] The intuition behind 
the algorithm is as follows. G p gets much smaller than G (some size T = n/Z % for a parameter Z), and 
there are at most Z % subinstances of that size. For each i, the subinstances correspond to a partition of P 
into consecutive pieces of size roughly n/Z l . Although each subinstance computes some detours which are 
short (only on Z nodes), since the subinstance has weights which are large (0(MZ 1 )), these short detours 
correspond to paths in G of large weight. 

After the computation of the short detours, the subinstance splits itself into a finer partition while in- 
creasing the maximum edge weight slightly in each new subinstance, and recurses on all the new subinstances. 
We show that the time due to the recursive calls is (^(ZMn^), provided Z = fl(\og n). The number of 
circumventing paths which are computed as candidates for replacement paths is at most 0(n 2 ) and so the 
overall runtime is 0{Mn u) ). 

Reminder of Theorem [3l The replacement paths problem for nodes s and t in an n-node graph with 
weights in {— M, . . . , M} can be solved by a randomized algorithm in 0(Mre u ) time if uj > 2, or in 0(Mn 2+e ) 
for every e > if uj = 2 . 

Proof. First compute the shortest path P = {s = V\ — > . ..,Ug = t} between s and t. Proceed as in 
Algorithm[TJ Step 3 of Algorithm Q] takes 0(n 2 ) time. We will later show that \L\ — 0(n 2 ). Because of this, 
the rest of the steps 6-16 of Algorithm Q] take 0(n 2 logn) time. 

Now consider Procedure^ Consider level i — 1, . . . , log n/ log Z of the recursion tree. It contains roughly 
Z l instances of size 0((n /Z z ) log 1 n). The weight size in any such instance is 0(MZ % ) For each of these, 
step 12 of Algorithm [2] takes 0((MZ l )((n/Z l ) log 1 n) u ) = 6((Mn u /Z iw - 1 ) log'" n) time, and step 15 takes 
6((MZ i+1 )((n/Z i )log i n) w ) = OdZMn^/Z^- 1 ) log*" n) time. 

Steps 16-19 add 0(n 2 /Z 21 ) detour distances to L for each of the Z l calls at the ith recusion tree level. 
The overall number of new items added to L at level i is 0(n 2 /Z l ). 

Consider Algorithm |3 For a call in the jth level of the recursion tree, and each subpath P!, steps 8-11 
take 6{{MZi +1 ){{n/Z3 +1 ) log J ' +1 n) w ) time via Yuster and Zwick's algorithm for SSSP [24]. The number of 
calls to steps 8-11 at the jth level of the recursion tree is Z J+1 so the runtime is 

M (n" log (j+1)w n)/Z^ +1 ^- 2 \ 

The number of new detours added to L in steps 8-15 is 0(n) for every level of the recursion tree. The 
number of levels of the recursion tree is 0(logn/logZ). The overall number of items in L at the end of all 
recursive calls is hence asymptotically 



log n/ log Z 

,logn/logZ+ V (n 2 /Z l ). 



This is 0(n 2 ) for any Z > 2, say. 

The overall running time is asymptotically 



log „yio g z ZMnU iu MnU (i+1) 
n 2 logn+ 2^ 



^i(w-l) 2(t+l)(u>-2) ' 

8=0 

If a; > 2, the above runtime is 0(ZMn u ), say for any Z > 21og"'^~ ' n. Hence picking Z = 
0(log CJ '^ CJ_ ' n) we obtain a runtime of 0{Mn u ). If u = 2, then we can pick Z — n e for any e > 0. 
The runtime becomes 0(Mn 2+2e ). 

We need to make sure that all random samples B work with high probability. The number of random 
samples B which are picked is at most 0(n log n) and the probability that a fixed one of them doesn't work 
is 1/ poly(n), so by a union bound, the probability that all of them work is at least 1 — 1/ poly(n). □ 



Algorithm 1 ReplacementPaths 



Input: G = (V,E), s, t; Global: L 

Compute the shortest path P = s = v\ —¥ V2 —>..., vi = t using Yuster-Zwick's algorithm 

Compute yvi,Vj G P their distance dp{vi,Vj) in P via dynamic programming 

Create a graph G' by removing edges on P and splitting every node v <E P into {w m , -y "*} 

Call RECURSIVE(G', P, s, i, M) 

Sort all elements (u,-, Vk, d) € L in nondecreasing order of d 

Store every (t>i,i) in a binary search tree T 

while I ^ do 

Pop (uj ,Vk,d) from L 

Find successor v q of Uj in T 

for every (v p ,p) G T with i < p < j do 
Set dfvpjWp+i] = d 
Remove w p from T 

end for 
end while 
Output d[-, ■] 



Procedure 2 RECURSIVE 



1 

2 

3 

•1 

5 

6 

7 

8 

9 

10: 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 



Input: (G',P',u,v,N), Global: L,P,d P ( v ) 

if G" has size O(Zlogn) then 

Compute the distances dc(-, •) between all pairs of nodes in P' in G' 
for every Uj, Wj e P', i < j do 

Compute of = dp(s, Vi) + dQ>{vi,Vj) + dp(vj,t) 
Add (vi , Uj , d) to L 
Return 
end for 
end if 

Compute matrix D so that (D * P)[u, u] is the distance of u and t> in G' 

Sample a set B of 0((|G'| Iogn)/Z) nodes from G' \ P' 

Replace all entries D[u,v] > N ■ Z with P[u, u] = 00. 

Compute D -k D 

for all Wj, Uj € P' with i < j and (P * D)[vi, Vj] < 00 do 

Compute d = dp(s,Vi) + (D -k D)[vi, Vj] + dp(vj,t) 

Add (vi , Vj ,d) to L 
end for 

Call ProcessSubpath(G', P', u, u, P, N) 
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Procedure 3 ProcessSubpath 



9 
10 

11 
12 
13 
11 
15 
16 
17 
18 
19 
20 



Input: (G', P', u,v, B, N), Global: L,P,d P (-,-) 

Split P' into Z consecutive buckets {P{, ..., P' z } of size \P'\/Z 
Remove all nodes in G" \ {B U P'} from G" 
for each P[ = ■ul i —>...—> t^ do 

Let P r and P' consist of the nodes to the right and left of P[, respectively 

Replace P r and P l by paths on \P'\/Z nodes as in Figure [3] creating graph G[ 

Find the distance di(u, v) between u and v in G\ 

Add (u' i ,v' i ,dp{s,u) + di(u, v) + dp(v,t)) to L 

Find the distances di(u,Vj) from u to all nodes Vj € P[ in G£ 

Find the distances di(vj,v) from all nodes Vj € P/ in GJ to v 

for all Uj e P[ do 

Compute d u = dp(s, u) + di(u, Vj) + dp(vj,t) 
Add (u^Vj^du) to L 

Compute d v — dp(s, Vj) + di(vj, v) + dp(v, t) 
Add (vjjv'^du) to L 
end for 

Remove all nodes corresponding to P' \ P[ from G' t 
CaU RECURSIVE^, P^u'^v'^NZ) 
end for 
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